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r-~^ . Abstract 

^^ . Consideration is given to a non-convex variational model for a shear experiment in the 

framework of single-crystal linearised plasticity with infinite cross-hardening. The rectangular 

Q, ] shear sample is clamped at each end, and is subjected to a prescribed horizontal or diagonal 

shear, modelled by an appropriate hard Dirichlet condition. We ask: how much energy is 
required to impose such a shear, and how does it depend on the aspect ratio? Assuming that 

rG , just two slip systems are active, we show that there is a critical aspect ratio, above which the 

jrt ' energy is strictly positive, and below which it is zero. Furthermore, in the respective regimes 

determined by the aspect ratio, we prove energy scaling bounds, expressed in terms of the 
amount of prescribed shear. 

> ! 1 Introduction and the main model 

r*~ I In this article, we study optimal energy scaling and microstructure formation in the elasto-plastic 

deformation of single crystals. The main model considered here is a geometrically linear, continuum- 
mechanical elasto-plastic energy based on the assumption of strong cross-hardening, i.e. that only 

^T ■ one slip system can be active at any given point of the material. The basic motivation for our 

investigation is the following question: can one devise an experiment to determine whether this 
single-slip condition, taken together with a surface energy which penalises geometrically necessary 
dislocations (GNDs), is a relevant constraint which needs to be factored into macroscopic models? 
The phenomenon of cross- (or latent-) hardening [U [^ [3] has the effect that activity in one 
slip system suppresses activity in all other slip systems, which leads to a loss of convexity in the 



JH ' minimization problem associated with the plastic deformation of single crystals [1]. In fatigue 

- -' experiments, as well as experiments involving only a single-pass deformation, a lamination behaviour 

is frequently observed [SJ HI [3 [5] , and this is believed to stem from cross- hardening [HIH]. Hence, 

similarly as in [5] , the interest in including cross- hardening in our model, as a side condition on the 

strain effected by plastic deformation. 

The surface energy we include corresponds to the core energy of GNDs. In the geometrically 
nonlinear setting, the density of the core energy would be of the form \{det F^)^^ {curl FP){F'p)'^ \ 
(see [nmnmi] for a discussion), where FP is the plastic deformation strain, but since we are working 
in the simpler geometrically linear setting, this leaves us with just the expression | curli^Pj instead. 
Our main results can be described heuristically as follows. We show that the inclusion of both 
cross-hardening and surface energy significantly affects the energetic scaling in a particular simple 
shear experiment for single crystals. The experiment modelled is that of a single crystal with a 
defined orientation of slip systems — realisable, for example, in a B2 crystal structure — and boundary 



conditions such that the cuboid crystal is rigidly fixed on its (square) top and bottom faces, according 
to a prescribed horizontal shear, while remaining free on its four vertical sides. Without either cross 
hardening or surface energy, the energy infimum in our system is zero (or equal to a certain amount 
of plastic work if a strain-hardening energy is included) when the aspect ratio of the crystal (base- 
square side-length/height) is below one, while the energy increases quadratically with the imposed 
shear magnitude when the aspect ratio is above one. This behaviour changes dramatically when 
multiple slip is forbidden and GNDs penalised: now the energy only vanishes for aspect ratios 
smaller than one half, while, again, for aspect ratios above one, the energy grows quadratically with 
the strain imposed by the boundary conditions. For intermediate aspect ratios, a new regime of 
linearly growing energy arises. Thus, we conclude that the experiment proposed here can be used 
to discriminate between those models with surface energy and cross hardening, and those without, 
and hence to determine whether the inclusion of these effects in macroscopic models for single-pass 
plastic deformation is physically reasonable. 

Kinematically, we restrict ourselves to the case of a cubic crystal structure with {011} (100) 
slip systems. This means that plastic deformation occurs only on planes with normal parallel to 
one of the three cube axes, and in the direction of one of two Burgers vectors lying diagonally in 
these planes. Furthermore, we include a surface-energy term consisting of the distributional curl 
of the plastic deformation, as well as a constant hardening. The domain is a cuboid given by 
rii = (0, 1) X (0, L) X (0, 1) C M'^, and we consider the problem of minimizing the energy 

El{u,P)^ [ \{Vu^(3%y,,\^ dx + a f \cml(3\+T f \I3\ dx (1) 

among vector- valued displacements u : O^ — > MP and matrix- valued plastic deformations /3 : O^ -^ 
]^3x3^ for non-negative coefficients a and r. Here, the subscript 'sym' denotes taking the symmetric 
part of the matrix in parenthesis. The first term in ([T]) is the linearized elastic energy of the specimen, 
and the second term penalizes GNDs, which are quantified by the row-wise curl of the Nye tensor. 
The third term represents a constant hardening, which for the most part will be neglected by taking 
T = (see, however, remark after Thm l3.l| ). 

The admissible displacements for the minimization problem are u G H^{flL), such that u{-, •, 0) = 
0, u{-,L, ■) = 7(1,0,0) for some parameter 7 > Qj. This models an experiment in which a cuboid 
specimen with clamped boundary conditions on the top and bottom undergoes a simple shear de- 
formation. Note that the free boundary conditions on the four sides of the specimen are essential 
in order to render this a physically realisable experiment, and that they also account for most of 
the technical difficulties in the analysis. While not subject to any boundary conditions, the plastic 
deformation /3 must be square integrable, and the absolute value of its distributional row-wise curl 
should be the density of a finite measure. In the sequel, lower bounds for the energy will be the 
focus of attention, and all test functions u and /3 will be (at least piecewise) smooth — this does not 
represent a further restriction, since any admissible deformation can be approximated in energy by 
a smooth one. 

The single-slip condition translates into a side condition on /?, which, in our proposed experiment, 
reads as follows. The single-crystal is oriented such that one of the cubic crystal axes points in the 
X3 direction of the specimen, and the other two lie in the xi — X2 plane, rotated by 45° with respect 



*For comparison, we also consider the diagonal- shear case, where the upper end of the specimen is fixed according 
to «(-,L, 0=7(1, 1,0). 



to the xi — X2 co-ordinate axes. Thus, 




ij ,s(x)( -l/^/2j(^|ll, (2) 

1\ /l/^/2\ 

-1 ,.s(x) I/V2 




for almost every a; = (a;i,X2,a;3) g fi^ and some coefficient s: Ol — > M. Note that there are indeed 
materials which exhibit such (primary) slip systems (or rotated equivalentijj) , in particular materials 
with B2 (caesium-chloride) structure, such as the intermetallic compounds Yttrium-Zinc [13j or 
Nickel- Aluminium. 

The article is organized as follows. In the next section we introduce a two-dimensional (sliced) 
version of our problem, and prove the key energy inequalities. Then, in section [3l we show that the 
relaxed three-dimensional problem can be reduced to the two-dimensional one, simply by adding 
up slices. In section 01 we show that a further reduction to a two-dimensional scalar model of 
plasticity does not, in our set-up with free lateral boundary conditions, yield reasonable energy 
bounds. Finally, some conclusions and open problems are discussed in section [S] 

2 The two-dimensional model 

We now introduce a 2-d model which basically corresponds to a vertical slice of the full 3-d model. 
Once the dimension reduction has been made, we proceed to obtain the desired upper and lower 
energy bounds. 

2.1 Description of the 2-d model 

Consider a rectangular specimen which occupies the domain Q.^ = (0, 1) x (0, L) C M^, and look for 
a vector displacement u: Q.^ — ?► M^ and a plastic-distortion tensor j3: VLl ^^ K.^^^ which minimise 
the linearised-plasticity functional 

El{u,P)^ [_ |(Vu-/3),y„,pda; + a / |V x /3|, (3) 

subject to the single-slip side condition 

m « ( J :J ) or (^ _J _J ) , (4) 

for a.e. x = (a;i,X2) G ^l- 

The curl term appearing in ([3]) models the energy of GNDs (with a interpreted as the line 
energy) , and is defined by analogy with the measure-theoretic quantity introduced in |12j , such that 



t There may also exist further Burgers vectors lying within the respective slip planes — wc show in section [3] that 
such considerations are irrelevant to our analysis. 



it is continuous with respect to a mollification of /3 which respects the side condition ^ . The most 
convenient way of explicitly making the definition is in terms of co-ordinates (^, 77) which are rotated 
at 45° with respect to the Xi, and hence aligned with the slip directions. Thus, for (3 e BV, say, 

|V X /3| = sup /_ 0,5^5,, + V',,,/3,,5 d^dr], (5) 

JnL 

where the supremum is taken over all 0, V' G C'q (fii, [—1,1]), and we have used the comma notation 
for partial derivatives. 

Note that, in such rotated co-ordinates, the side condition (j4]) takes the simpler form 

/3(?) '?) « ( 1 n ) o'" ( n n ) ' (6) 



^ 1 y V ^ 

for a.e. {^,i]) £ ^l , which will be extremely useful in the sequel. Also note that, for smooth /3, ([S]) 
is just the L^-norm of the ordinary row-wise curl of /3. 

Our main goal in what follows will be to obtain lower bounds on E^ as the reciprocal aspect 
ratio, L, is varied, subject to two sets of boundary conditions, corresponding to diagonal and hori- 
zontal shear, respectively. The infimum of El will be taken over (u, /?) which satisfy the boundary 
conditions under consideration, and we will write Jl ~ inf i?L(', •)■ 

The boundary conditions to be treated are, when expressed in Xi-co-ordinates, 

(BCl) u = 7(1, 1) at X2 = i and u = (0, 0) at X2 = 0, 

(BC2) u = 7(1, 0) at X2 =: L and u = (0, 0) at x^ == 0, 

where 7 > is a measure of the fixed average (diagonal or horizontal) shear. 
We begin with the following pair of elementary results: 

Proposition 2.1. Under BCl or BC2 , Jl is a monotonically decreasing function of L. 

Proof. Fix L > 0, and choose < 5 < L. Next note that we can write flL-s — ("^x^) ^l U A, 
where meas(A) > 0. Also note that El is invariant under an overall scaling of fi^, via u{Xx) = 
u{x),${Xx) = jP{x), for a scale-factor A, say. The result follows by the super-additivity of the 
energy as a set function. D 

Proposition 2.2. Under BCl or BC2 , Jl is a left- continuous function of L. 

Proof. Fix e,i > 0, and take a smooth pair (w, /3) on Q,l such that El{u,(3) < Jl + e. Now scale 
u and (3 as in the above proposition to get (m,/3) defined on 5 = (t^) ^l C i^L-s- Extend (u,/3) 
to the whole of ^l-s by constancy in the a;i-direction, thus preserving the boundary conditions at 
a;2 = and X2 = L ~ 5 . This construction induces no extra curl on xi = 1 — -j-, and the energy in 
O.L-5 \ S tends to zero as (5 — > 0. Hence lim5_j.o+ Jl-s < Jl + ^- Since e was arbitrary, we get the 
desired result. D 

2.2 Diagonal shear 

While a diagonal shear would be rather difficult to impose in a real experiment, due to slippage, 
it turns out that most of the test-function constructions required for horizontal shear have simpler 
analogues in the BCl case which are much easier to explain in the first instance. For this reason, 
we now present the BCl model by way of an introduction to the more complex BC2 problem. 




Figure 1: A minimising pair for BCl in the case L > 1. 

In terms of rotated co-ordinates, {£,,'']), aligned as in Figure [U the Dirichlet condition of BCl 
at the top boundary becomes (u^,u,,) = 7(0, v2), while, of course, (m^jU,,) = (0,0) at the bottom. 
Using this, it is easy to construct zero-energy test functions whenever L > 1. 

Theorem 2.3. Subject to BCl, L > 1 ^ ,/l = 0. 

Proof. Take a simple-shear band S, of width e, which just misses the horizontal boundaries, as in 



Figure [TJ Let u^ 



P ~ below S, u^ = 0, u,-i ~ \/2^, (3 — above 5", and interpolate linearly 



between the two to ensure that 



, 



(7) 



on S', in (^, 77) co-ordinates. Clearly the boundary conditions are satisfied by the resulting test 
functions, and it is trivial to check that El{u,(3) =0. D 



Theorem 2.4. Subject to BCl, L < 1 ^ 3 cl > such that 

^(1 -L)<Jl< min }^,J-{1-L)+ clg^ 



(8) 



Proof. For the lower bound, fix L < 1, and choose a pair (m,/3) satisfying BCl. By ([3]) and ([5]), we 
immediately get 



drj 



L^(nL) 






<El{u,(3). 



(9) 



L^(nL) 




1 -L 



Figure 2: The line l^ and the region R used in the lower-bound proof of Theorem 



Referring to Figure [2l we fix ^ such that the dotted line l^ (of constant ^) lies in the region 
R C r^L which is bounded by the two solid diagonals. Now integrate -^ along l^ from u^, = to 
u^ = \/2j to get 



and integrate over all such ?^ C i?, to arrive at 

Using \R\ = L{1 — L), along with the usual L^ ^^ L^ embedding, we get 



7(1 -L)< 



dUr. 



dr] 



<{L(l-L))^ 



LHR) 



dUn 



drj 



(10) 



(11) 



(12) 



L^[R) 



Hence, by ©, El{u,/3) > 7^(1 — L)/L, as required. 

For the upper bound, first note that the first member of the upper-bound set is attained by the 
minimising (linear) purely elastic deformation. Next, we get arbitrarily close to the second member 
with the following construction, which is depicted in Figured 

Divide ^II into five regions: Re (a slightly narrower version of R from Figure [2]), two transition 
layers of width e centred on the corner-diagonals of constant ^, and two triangles which surround 
the i?e /transition- layer sandwich. 

Set u^ = on fi^, and, on R^, determine m^, for each fixed ^, by linearly interpolating between 

the Dirichlet conditions at X2 = and X2 = L, and set /3^j = -qt-- On the upper and lower triangles 

take /3 = and, respectively, Mj, = \/27 and u,-, = 0. This gives a total elastic energy away from the 

2 . — 

transition layers amounting to ^(1 — L — v2e), and zero plastic energy. 




Figure 3: Construction used in the upper-bound proof of Theorem 



On the transition layers, {u, (3) can be chosen such that the elastic energy is small, and the plastic 
energy is 0{aj). To be specific, divide the upper transition layer into a small triangle, A, in the 
bottom left corner of Ol, and a narrow strip, B (see Figure |4^). Then, by analogy with the proof 
of Theorem 12.51 (see below), we can make the energy in A arbitrarily small, as follows. Thus, take 
a function m£($) with (piecewise) ^"-shape, representing a steep transition from to ^ along the 

diagonal boundary of A. Then set (3 = Vm^, put v = u,, — u^ and minimise ||Vw||l2(^), subject to 
w^ = on a;2 = and m^ = m(^ on the diagonal boundary of A. This can be done by reflecting A 
about its vertical boundary and solving the Laplace equation for v on the resulting domain, which 
gives arbitrarily small energy as a — > 0. 

Meanwhile, on B, interpolate u,, between the u^-profile at the boundary with A and the Dirichlet 

condition at X2 = L, for each fixed ^, and accommodate the resulting twist by setting /3^^ = -of-- 

The result of this is that we have, for small enough a(e), 0{€) elastic energy and 0(0-7) plastic 
energy on the upper transition layer. Clearly, an analogous construction works on the lower transition 
layer, and hence, letting e — > 0, the desired upper bound follows. 

n 

An interesting question is whether J^ jumps from zero to a strictly positive number as L decreases 
through unity: we now answer this question in the negative. 

Theorem 2.5. Subject to BCl, Ji = 0, and hence, by Provosition \2. 2[ Jl is continuous at L — 1. 

To prove Thcorcm l2.5[ we need the aid of the following Lemma. 

Lemma 2.6. The family of functions fa{x) = 1 — x", a G [0, 1], lies in iJ2(0, 1), and \\fa\\ 1 > Nj 

as a — > 0+. In particular, if 2 (0,1) is not embedded in L°°(0, 1), since WfaWL-^io,!] = 1 for 
a G [0,1]. Moreover, the reflected functions fa{x) — 1 — {sgn(x)x)°' lie in H^ (^^1,1), and also 




u,(^) 



Figure 4: Construction used in the upper-bound proof of Theorem 12.41 (a) The upper transition 
layer, comprising regions A and B, (b) Profile of itr, on the diagonal boundary of A. 



Proof. Clearly, /„ e L^, ||/q||l2 \ and ||/a||L~ ~ 1- Now, turning to the H2 semi-norm, we 
define 



for ia,x,y) € (0,1)3. 

Thus, Q(0) = 0,Q(1) = 1, and P(0) = 0. Moreover, 



P(a,x,2/) = (x-"-2/")^ 



(13) 



dP 
da 



2{x" - y°'){x°' Inx ~ y'^lny) > 0, 



(14) 



so that P{-,x,y) is increasing, and clearly P{a,x,y) — !• P{0,x,y) = as a ^ 0+. Thus, by 
dominated convergence, (5(a) — ^ as a — >■ 0"*", as required. 

For fa, the semi- norm integral can be split into two 'diagonal' terms, which can be handled as 
above, and two 'cross' terms, given by 



/.I 



1 Jo 



(x" - (sgn(y)y)")^ 
{x-yy- 



dxdy = 



1 (-l /^Q -,a\'2 



"'0 



(x" - z")2 
{x + zf 



dxdz 



(15) 



where we changed variables via z = —y in the final equality. Clearly, the last integral is smaller than 
(5(a), and increasing w.r.t. a, which gives the required result. D 

Proof of Theorem \2.5i Here it will be convenient to align the rotated (^, rj) co-ordinates as in Figure 
[SI with the origin taken to be at the centre of fii. 

Given this, consider the sigmoid function ga{£.) defined by 



ffa(0 = ^(l-sgn(0(sgn(0V2e 



(16) 



on the interval (— 1/V2, l/-\/2), which represents a steep transition from to V^j (see Figure [6]). 
Now define, in the (^, 77)-frame, u^ = (0,5q(^)), and put 



/3 = Vu' 



P 







(17) 




Figure 5: Orientation of the {£,, if) co-ordinates on fii used in the proof of Theorem 




Figure 6: The function ga{C} used in the proof of Theorem 12.5 




Figure 7: Shape of the Dirichlet data on fJi used in the proof of Theorem 12.5 



Thus, /? is curl-free, and, for a general test function u, we define v = u — u^ and set v^ = 
(44> u^ ~ 0, consistent with BCl). The aim is now to show that the minimum of ||Vw,J|i2, subject 
to the boundary conditions 



Vni^,v) 



V2j- 






X2=0 
X2 = 1 



(18) 



tends to zero as a — )- 0: this will be achieved by using the continuity of the (Dirichlet) solution 
operator to Laplace's equation, Avri = 0. 

In order to do this, it is convenient to consider reflected (about xi ~ 0) Dirichlet data on the 
double of ill, with resulting lateral boundaries then identified to make a cylinder. Thus, we work 
on the manifold ili = S^ x (0, 1), with boundary §^ U S^, and the Dirichlet data for Vrj on each 



boundary §^ have the same shape as fa from Lemma 12.61 (see Figure [7]) . By the symmetry of this 
data, and invariance of Aw^ = under reflections, we see that solving the Dirichlet problem on f2i, 
for a function still denoted by w,,, gives a solution of the mixed Dirichlet/Ncumann problem on f2i. 
Now, Proposition 1.7 on p. 307 of [13] tells us that the solution map to the Dirichlet problem on 
0,1 takes H^i^dOi) continuously into iJ^(rii), while our Lemma [2.61 gives, modulo trivial scalings 
applied to /„, ||Tr(w„)|| i - -> as a — ?► 0, as desired. D 

Remark 2.7. In the above construction, if we use a piecewise-linear transition function with 3-pieces, 
rather than g^, then we get an 0(7^) amount of energy, as can easily be checked, using the trace 
inequality. 



2.3 Horizontal shear 

In this section we minimise El subject to the horizontal-shear condition BC2. In a (^,77)-frame 
aligned as in Figure [U the Dirichlet condition dX X2 = L is {u^,Urj) = 7(— 1/a/2, l/-\/2), while 
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Figure 8: Zcro-cncrgy construction for Theoreni l2.8 



(w^, u^) = (0, 0) at X2 = 0. 

Theorem 2.8. Subject to BC2, L>2^ Jl=0. 

Proof. Wc use basically the same construction as in the proof of Theorem 12.31 except that here, 
since L > 2, there is enough room to fit two alternating, non-intersecting shear bands into il^, each 
of which misses the Dirichlet boundary (see Figure [8|) . 

To be explicit, we take two shear bands, given by Si — {xi +X2 < L} O {xi +X2 > 1 + L/2} and 
'S'2 = {3:2 — xi > 0} n {x2 — xi < — 1 + L/2}, which just touch at the lateral boundary xi = 1. Above 
Si set u^ = — 7/-\/2, u,, = 7/\/2, /? = 0, between 5*1 and ^2 set u^ = 0,Urj = 7/V2, /5 = 0, and below 
^2 put U(^ ^ Uri ~ 13 = 0. Then linearly interpolating u across Si and S'2, and requiring that Vm = j3 
on 5*1 U S'2 clearly gives a pair of test functions which satisfy BC2, and for which El =0. D 

A certain amount of energy is required to produce a horizontal shear when L is sufficiently small, 
which is the content of the following two theorems. 



Theorem 2.9. Subject to BC2, L <l 

„2 



3 cl > such that 



r.^l-L)<JL<niui{^,^il-L) + CL., 



(19) 



Proof. The proof is very similar to that of Theorem 

Thus, 
[21) , to get 



Thus, for the lower bound, take a test pair (u, /3), and integrate -^ over the strip R (see Figure 



dUrj 

dr] 



> 



L2 



4L 



(20) 
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Figure 9: Transition layers used in the proof of Tlim l2.9 



and integrate -g/ over an analogous strip pointing in the ^-direction to get the same estimate for 
-^ . Hence, El{u,I3) > 7^(1 — L)/2L, as required. 

For the upper bound, it is easy to see that the first member of the upper-bound set in (J19p is 
attained by the minimising purely elastic deformation, while the second member can be approached 
by a construction similar to, but more complicated than, that used in the proof of Thm 12. 4[ the 
difference being that here we need two transition layers for each of m^ and m^, and we have to take 
a little care on the regions where they intersect. 

Thus, referring to FigurelHl "^c define, for i = 1, 2, the transition layers Tj' to be narrow strips of 
width e, centred on the corner-diagonals of constant ^, and Tf to be strips of the same width, but 
centred on the corner-diagonals of constant 77. Between the r/*, u^, is obtained by linear intepolation 
between the Dirichlet conditions, in the 77 direction, while above and below the T^'^-sandwich, u^ is 
set equal to the appropriate Dirichlet condition. The other component of the displacement, mj, is 
obtained on the complement of the Tf in an analogous manner. 

Next, Uri is defined on T^ as follows. Divide the part of T^ below the intersection with T^ into 
two parts, A and B' , such that A is the same as in Figure [3] and B' is a truncated version of B (see 
Figure ITU)) . Define m,, on A in exactly the same way as in Thm 12.4] with a > small and fixed, and, 
for each fixed ^, interpolate from the diagonal boundary of A up to the unique linear profile on the 
upper constant-?7 boundary of B' which connects u^^ continuously across T]^ . On the remainder of 
T]' , define u,, by interpolating up to the Dirichlet condition at X2 ~ L, for each fixed ^. 

On T2 , Uri is determined by a trivial symmetry transformation applied to m^ on T^ , and likewise 
for w^ on the Tf . For the plastic distortion, let /3 = Vu^ on A, set (3 = Q onT^ f)Tf , and on the 
rest of T|' let /3^j — -^ . Define /? on the remaining transition layers by symmetry, and otherwise 
by appealing to Figure [9l 

As a result of this, there is very little elastic energy on the transition layers, while the plastic 
energy is concentrated there (particularly on the boundaries). More precisely, we have 11-^11^2 = 



12 




Figure 10: Subsets of T"^ used in the proof of ThmEJ 



0(e) on T^ and ||-^||^2 = 0{e) on Tf , provided a(e) was chosen small enough, while, by the 
symmetry of the construction, -g^ + -^ = ^(l)i pointwise, on the transition-layer intersections. 



Also note that, by symmetry, -^ + -^ = on the central square, where /3 = 0, and thus there is 
no off-diagonal elastic energy there. 

Finally, it is easy to see that all the contributions to the plastic energy are 0{a^), and hence, 
for our test functions, 

„2 



EL{u,P) = ^{l-L) + 0{a^) + 0{e), 



so that, letting e ^^ 0, the desired upper bound follows. 



(21) 

n 



We have to work a little harder to prove strict positivity of the energy in the intermediate case 
1 <i < 2. 



Theorem 2.10. Subject to BC2, 1<L<2^3cl>0 such that 

2L 



clcjT 



O + A/cr2 ^ 2CL72 



< Jl < min{ -^^,2\2'^a 



(22) 



Proof. For the lower bound, we combine a recent generalisation of Korn's inequality from |15| which 
is nicely adapted to our energy functional with a Lemma from |12j which allows us to control (one 
component of) j3 via the curl, provided the single-slip condition is satisfied. 

First, let 1 < L < 2, and, for test functions (u,/3) satisfying BC2, we define the average- 
symmetric tensor A via 



^(u,/3) = (Vw-/3)- 



ni 



(Vw - P) 



skew 



dx. 



(23) 
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Figure 11: The region Q used in the proof of Theorem 12.101 



Then, by Theorem 11 of fT51, there exists a constant C{^l) such that|3 

/ |A,kcw|' dx<c(f_ |^y,n|' dx + II V X ^|!il(o^) j , 

and hence there exists a constant, skew-symmetric matrix K{u,l3) such that 

\iVu -13)- Kiu,(5)f < CEl{u,P) (l + :^£%^ 



(24) 



(25) 



Now, consider the square Q C f^L, of side 4=, which is bounded by the hues {xi + X2 = L/2}, 
{x2 — xi = — 1 + L/2}, {xi + X2 = 1 + L/2} and {x2 — xi = L/2} (See Figure [TT]): for convenience, 
we have put the origin of the {£_, 77)-co-ordinates at the bottom corner of Q. 

By analogy with the proof of Lemma 4.3 of [H], we define the fohowing subsets of Q: 



ui 



(?) 
Xv) - 



{e e (0, 1) : /3^4(e, v) ^ for a.c. ij e (0, 1)} , 
= {?; e (0, 1) : /34^(e, v) ^ for a.c. f £ (0, 1)} , 



(26) 
(27) 



and, from the side condition Q, we conclude that at least one of w'^^ and a;''') must be a null set. 
For the remainder of the proof, we may assume w.l.o.g. that a;''" is a null set. 
Let 

P = nin||<xi+X2<l + ||, (28) 

so that P is, so to speak, the extension of Q in the ^-direction. 



-fNote that, at the time of writing, l|24|l is only known to hold in 2-d, although analogous inequalities in higher 
dimensions have been conjectured by Neff et aL|16|. 
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C- 



Then, since uj^^' is null, wc obtain 
and it follows from ©, (HSI, ^^ and L^ ^^ L^ embedding that, 



du/L 



di 



<CEl 



LHP) 



and 



dri 



K, 



^v 



< CEr 1 



L\P) 



< CEr 1 



El 

El 
a2 



If^inWLi 



LHP) 



Define strips Ti, Ta C P by 



Ti^Pn{xi+X2<l}, T2 = Pn{xi+X2>L}, 



(29) 

(30) 

(31) 
(32) 

(33) 



such that Ti intersects the bottom boundary, and T2 the top boundary, in a set of length 1 — L/2. 
The strategy now is to show that 



Wu^hnn) < CEl, 

|w5 + 7/V2|Ui(T,) < CEl, 



(34) 
(35) 
(36) 



Then, integrating -^ between Ti and T2 will give a lower bound on K^,j, and a contradiction if El 
is too small relative to the average shear 7. 

We begin by integrating -gr- from (Co(^)j'7) on the a;2 = boundary of Ti to the point (^,77) in 



the interior (see Figure [T^ . to get 

which, upon integration over {^,1]) G Ti, gives 



// \u^\dCdr, < ff f 



iv) 



du^ 



d^ 



ie,v)ded(dTj 



< C 

< CE 



du^ 



dc 



LiiQL) 



L' 



(37) 

(38) 
(39) 
(40) 



and therefore ([M)) holds. The estimate ([55]) follows by an analogous integration on T2. 

For the upper bound on K^,j, consider the small triangle i?i = Ti n {xi — X2 > ^} next to the 
a;2 = boundary (Figure [T^b)- 
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(^.11) 



(^o(Tl).Tl) 

Figure 12: Construction used in the proof of Eq. (p^ 



^.11 2©) 



(^.Tll©) 




(b) 



©(ll).Tl) 




(^.Tlo) (4l(Tl)>Tl) 



Figure 13: Constructions for the lower (a) and upper (b) bounds on i^^^ in Thm [2?T0l 
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Integrate -^ from X2 = to the point (^, 77), holding ^ fixed: 



ud^,v)=f^itri')dv', (41) 

then integrate (|1T|) between ^1(77) and £,2(11)^ 3'iid sum over aU Une segments (^i,'C2) lying in ^1, to 
get 

// u^c^v) d^dr^ ^ r r" r ^(^^^') ^^'^^^^^ (42) 

where [771,772] = n^(i?i), and 11,, denotes projection onto the 77-axis. 

Adding and subtracting K^ri to the integrand on the rhs of (|42p . and using ([32]) and (p4| . we 
arrive at 

\K,,\<CEi{l+(l+'^-^y\- (43) 



For the lower bound on A'^,,, we work on the square Q, and, in particular, the subsets Q n Tj 
thereof (see Figure [T3b). 

Firstly, we integrate -^ between the two points (^, 771 (^)) G Q n Ti and (^, 772(0) G Q n T2: 



/ ^(^,^')d^' = «^(^,^2)-««(0»?i). (44) 



^'72(0 Q^ 

Then, we integrate (|44|) over ^ e nj((5), 771 e n,,(ri) and 7^2 £ n,,(T2), to arrive at 



^(C, 7,') diid£,drjidin = -^ (^1 - I) (^1^ uj derf^2 - 1^ ^c dC*?i) • (45) 

Thus, adding and subtracting A'j,, from the integrand on the Ihs of pS)) . and using (|5^. ([M|) 
and (|35|) . gives us 

|A^«.I > C f^ - si fl + (^1 + ^y "l V (46) 

which, together with (|43p . implies 

4(l+(l + ^)')>C'7. (47) 

Squaring, and using (a + 6)^ < 2(a^ + 6^) plus the quadratic formula, leads to the required lower 
bound on J^. 

For the upper bound on J^ , we note that the first member of the upper-bound set in (|2.10p is 
attained by the minimising purely elastic deformation satisfying BC2, while the second member is 
attained by a simple crossing-shcar-band construction, as follows. 

Referring to Figure [T4l inscribe, in ^l, crossing shear bands 5*1 and ^2 which just miss the Dirich- 
let boundaries, and which have width -j^{L — 1). Then define uj and Urj by linearly interpolating 

between the boundary conditions across the respective Si. Finally, set /3 = Vu on 51^ \ {Si n ^2) 
and /3 = on S'l n S'2. An easy calculation shows that the elastic energy vanishes, while the plastic 
energy is given by 2v270', as required. D 
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(a) 



(b) 





Figure 14: Upper-bound construction for Thm 12.101 (a) in the reference configuration, and (b) the 
deformed configuration. 



Remark 2.11. If the single-slip condition is dropped then J^ = for L > 1, by essentially the same 
construction as in Figure [Ml the difference being that we may take /3 = Vm on the whole of ^l, and 
not just away from the shear-band intersection. 

Finally, we show that, as in the diagonal-shear case, the transition from zero to positive energy 
is continuous. 

Theorem 2.12. Subject to BC2, L = 2^ Jl = 0. 



Proof. The proof is similar to that of Theorem [ 

Take L = 2, and let the (^, 77)-co-ordinates be centred at the mid-point of the right-hand lateral 
boundary, as shown in Figure [T51 

We construct test functions (m^,/3) as follows. On {^ > 0,7^ < 0} let (u^,u^) = 7(0, l/\/2) and 
/3 = 0, on {5 > 0, 77 > 0} let 



7 
V2 



(^/2r,)^ <-^, /3 = V7.^ 



and on {^ < 0, 7? < 0} let 



«f=0, ^^ = ^ (1 - (-^/2Cr) , /3 = Vu' 



V2V 



(48) 



(49) 



for a e [0,1]. 

Clearly, V x /3 = 0, and, writing v = u — u^ a.s usual, we wish to make El = |l(Vw)sym|l/,2 
arbitrarily small by letting a — > 0, subject to the boundary conditions 



(^'?:«.)) 



0'^ 



^{{V2rjr-1),0] : X2 = 2, 
{-l + {-V2Cr)] : 3:2 = 



Since we have 



ll(V7;).y„,|li. 
it is enough to show that ||Tr(z; 
diagonal-shear case. 



<C{\\Vv^\\h 



Wv. 



nWL^) ; 



(50) 



(51) 



ffi(ao2) 



as a ^ 0. The proof of this follows just as in the 

D 
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L=2 




Figure 15: Orientation of (^, 77)-co-ordinates used in the proof of Theorem 12. 121 

3 Extension to 3-d: B2 crystals 



A sceptical reader might be tempted to dimiss our 2-d model as somewhat artificial, but in fact the 
results of section [5] carry over almost immediately to the case of crystals with B2 symmetry. 

As described in the introduction, two of the slip-plane normals in the B2 crystal are aligned with 
the (^, 77)-axes, as before, and the third is aligned with the C (resp. a;3)-axis, which is taken to be 
pointing out of the page. 

In co-ordinates Xi aligned with the shear-sample faces, the domain occupied by the crystal is 
taken to be fiL = (0, 1) x (0, L) x (0, 1), and we impose a horizontal shear of magnitude 7, such that 
the boundary conditions take the form 



7(1,0,0) at X2=L, w^: (0,0,0) at 2:2 = 0. 



(52) 



For the purposes of deriving the required lower energy bound, the single-slip side condition on 
the plastic distortion /3: fl^ -^ M.'^^^ is now taken to be, in the (^, 77, C)-frame, 



,33 




(53) 



for a.e. (Ci^iC) G ^II, and there is no restriction on the 3^(^,77, C). Note that this is a somewhat 
weaker condition than the one in ([2]) , in the sense that we allow any Burgers vector which lies within 
a given slip plane, and not just those in the direction of the cubic-latticc-face diagonals. Of course, 
any lower energy bound for the relaxed condition still holds for the original one. 
In this setting, for a.e. (^, ry, C), we have 



/3?r, /3«c^ 
/3 = I /3r,5 /3„c 



(54) 
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and for a.c. slice of constant C (call it S"^), the 2-d side condition holds: /3^j, = or /3,,^ = for a.e. 

Now, consider the energy functional ([T|), with r set to zero and subject to the relaxed side 
condition ([5^ . Denote by (Vw — /3)s™ ^^ic (^, ?7)-coniponents of the elastic strain, and, for fixed (, 
define the reduced slice energy, Ej^, by 

Eiiu, /?) = // (Vu - (3)1^^^ ' + \d^(3^,\ + \d,M d^dr,. (55) 



Then, clearly, El{u,(3) > L Ej^ dC,, and, since all of our 2-d lower bounds from previous sections 

apply to i?^, for a.e. C> they also go through for the B2 crystal oriented as above. Meanwhile, the 
upper bounds for the energy can be recovered in the following way. First, simply take the same test 
function, (m, (i), as for the 2-d case, extend constantly in the C, direction, and let /3f . = (3.^ = 0. This 
establishes the upper bounds subject to the relaxed side condition ([53)) . The original side condition 
can be recovered by introducing fine oscillations between the two possible Burgers vectors, in the 
active slip planes along the slip-plane normal. Note that such a microstructure is a gradient field 
even without the introduction of a lattice rotatioi^^, and thus does not make any contribution to the 
curl term. 

In summary, we have proved 

Theorem 3.1. The upper and lower bounds of Theorems \2. S'UKTB also hold for the energy (Qp with 
T = 0, subject to h52j) . as applied to a B2 crystal occupying the domain D,l — (0, 1) x (0, L) x (0, 1) with 
slip-plane normals oriented as above, and such that the cross-hardening condition i53\) is satisfied. 



Remark 3.2. For the case of a non-zero, constant strain hardening, i.e. r > 0, the energy scaling 
is as follows. For i > 2, the minimal energy is now at most T7, which is easy to see by inspecting 
the 2-d test function used previously. For L < 1, we still have the same quadratic scaling of the 
energy with respect to (large) 7 as before, since the upper-bound test function of Figure [S] gives a 
strain-hardening energy of order T7. In the regime 1 < L < 2, the linear scaling of the energy with 
respect to 7 is preserved: the proof for the lower bound goes through unchanged, while inspection 
of the test function for the upper bound yields an additional contribution to the energy of order T7. 
Furthermore, for a physically realistic model, we would expect r to be small compared to both the 
elastic modulus (implicitly equal to unity here) and a. This is due to the fact that the L^-pcnalty 
on 13 only accounts for self-hardening, i.e. the hardening which occurs if the specimen is deformed 
purely in single slip — all other hardening is taken care of by the single-slip side condition and the 
surface energy. 

4 The scalar model 

In this section, we take a brief (and salutary) look at an ostensibly simplified scalar version of the 
model ©-([S]) which is analogous to that treated in section 4 of [H]. Rather than imposing a soft 
Dirichlet condition on the whole of di^L as in [12], however, we continue to use a hard Dirichlet 
condition to model a fixed horizontal shear imposed by clamping both ends of the single-crystal 
specimen. 



/O /3i2(x2) 0\ 

§It is in fact a rotated version of in Zi-coordinates. Note that the row-wise curl of such a 

\0 0/ 

matrix vanishes identically, independently of the scalar function fli2{x2)- 
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L = 




u = 

Figure 16: Test functions used in the proof of Theorem 14. II 

As before, the rectangular specimen occupies the domain H.^ = (0,1) x (0,i) C M^, but this 
time we focus attention on just the xi-component of the deformation field, denoted by u: fl^ — > K, 
and the plastic strain, /3: Hl — > K^, which are to be obtained by minimising the functional 



Ei{n,P) = 



\Wu-l3\^dx+ / |V X /3|, 



(56) 



subject to the side condition /3i = ±/32, in the sense of distributions in Q^, and where the second 
term in (|56p is the measure-theoretic quantity defined in |12| . This quantity can be obtained from 
([5|) by making the substitutions 13^,] -^ l^ri and /3^j -^ /3^ in the usual rotated (^, 77) co-ordinates. In 
these co-ordinates, the side condition on /3 takes the form: /3^ = or /3^ = for a.e. {£,,t]) G Ul- 

We prescribe free boundary conditions on the two 'vertical' sides of f2i, and a fixed average 
shear, 7, in the xi direction on the remaining two sides. In other words, (n, Vw — /3) == on 
dflL n {xi = 0, 1}, where n is the unit normal to dflL, u = on dflL Ci {x2 ~ 0} and u = 7 on 
OUl n {X2 = L). 

Now we write J£ = inf i?£(-, •), and ask how J£ varies with the L. It turns out that analysing 
this question is no easier than in the full vector-valued case, and that, in fact, there is a qualitative 
difference between the results for the scalar model and those for the vector-valued model subject to 
either BCl or BC2. 

First of all we note that J£ is left-continuous and monotonically decreasing, by the same ar- 
guments used in Propositions 12.11 and 12.21 Then the next Theorem, combined with the foregoing 
remark, implies that J£ = for i > ^ and that there is no jump in the energy at L = i: the 
difference between this result and those for BCl and BC2 is explained by the fact that there is no 
constraint on the vertical deformation in the scalar model, which is rather unphysical. Whether the 
energy is strictly positive for L < -^ remains an open question. 

Theorem 4.1. Subject to u ^ at X2 ^ and u ^ j at X2 ~ \, we have Jf = 0. 



Proof. Let i = i and centre the (^, 77) co-ordinates at the mid-point of the top boundary, as shown 
in Figure \W[ 

Now, for each a € (0, 1), construct test functions {up, (3) by setting m^ = /3 = on {^ > 0, t] < 0}, 



Up 



7(2V2C)", /3 = Vw/j on {^>0,r/<0}. 



(57) 
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infEiii 



Quadratic 
\Scaling 




Zero Energy 



Figure 17: Schematic of the different scahng regimes for the energy ([T]) with r = 0. The sohd hne 
represents the energy when both tr > and the single-slip condition is enforced — in this case, the 
energy vanishes for L > 2, scales linearly with the (large) applied strain 7 for 1 < L < 2 , and scales 
quadratically for L < 1. The dashed line shows the energy when cither ct = or the single-slip side 
condition is dropped. In this case, the linear energy-scaling regime is no longer present. 



and 



Thus, V X /3 
deformation, and we look for v which minimises ||Vw||?2, subject to the Dirichlet conditions 



U0 = ji2V27j)" , l3 = Vup on {C<0,77>0}. (58) 

0, and, in the now familiar way, we write v = u — up, where u is the unknown 

L2j 





7(1 - (2V2C)") 
7(1 - (2V277)") 



X2 
X2 
X2 



r/<0 



(59) 



The fact that ||Vi'||i2 can be made arbitrarily small as a 
of Theorem 1 2. 5 [ given the shape of the Dirichlet data. 



follows by analogy with the proof 



n 



5 Conclusions 

We have derived energy-scaling results and explicit bounds for the elasto-plastic energy in a real- 
istic single-crystal shear experiment. The different scaling regimes under various assumptions are 
summarised in Figure 1171 As discussed in the introduction, the results obtained should provide 
a means of discriminating between two different classes of models in crystal plasticity: those with 
strong cross-hardening and surface energy, and those without. It is, however, open as to whether our 
results carry over from the B2 case to more common crystal structures, such as face-centred-cubic 
or hexagonal-close-packed lattices. 

Of course, in this article, we have only considered a geometrically linear elasto-plastic energy, and 
while our upper bounds should still hold in the geometrically nonlinear case, using a multiplicative 
decomposition of the strain, obtaining the lower energy bounds is a much more challenging problem. 



22 



and is the subject of ongoing work. A further objective is to carry out numerical simulations of the 
deformation predicted by our modelling. 

Finally, the experiments themselves are also work in progress, and will be performed in collab- 
oration with the Max Planck Institut fiir Eisenforschung in Diisseldorf — one of the main practical 
problems in this regard will be to impose the boundary conditions without slippage or other spurious 
deformation. 
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